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In the context of the Randall-Sundrum (RS) single-brane scenario, we discuss the bulk geometry 
and dynamics of a cosmological brane in terms of the local energy conservation law which exists 
for the bulk that allows slicing with a maximally symmetric 3-space. This conservation law enables 
us to define a local mass in the bulk. We show that there is a unique generalization of the dark 
radiation on the brane, which is given by the local mass. We find there also exists a conserved 
current associated with the Weyl tensor, and the corresponding local charge, which we call the Weyl 
charge, is given by the sum of the local mass and a certain linear combination of the components of 
the bulk energy-momentum tensor. This expression of the Weyl charge relates the local mass with 
the projected Weyl tensor, E^, which plays a central role in the geometrical formalism of the RS 
braneworld. On the brane, in particular, this gives a decomposition of the projected Weyl tensor 
into the local mass and the bulk energy-momentum tensor. Then, as an application of these results, 
we consider a null dust model for the bulk energy-momentum tensor and discuss the black hole 
formation in the bulk. We investigate the causal structure by identifying the locus of the apparent 
horizon and clarify possible brane trajectories in the bulk. We find that the brane stays always 
outside the black hole as long as it is expanding. We also find an upper bound on the value of the 
Hubble parameter in terms of the matter energy density on the brane, irrespective of the energy 
flux emitted from the brane. 

PACS numbers: 04.50.+h; 98.80.Cq 



I. INTRODUCTION 



The braneworld scenario has attracted much attention in recent years [1] . In this scenario, our universe is assumed to 
be on a (mem)brane embedded in a higher dimensional spacetime. There are many models of the braneworld scenario 
and corresponding cosmologies. One of them that has been extensively studied is the braneworld cosmology based 
on a model proposed by Randall and Sundrum [2] , in which a single positive tension brane exists in a 5-dimensional 
spacetime (called the bulk) with negative cosmological constant, the so-called RS2 model. In this paper, we focus our 
discussion on this single-brane model. 

In many cases, the 5-dimensional bulk geometry is assumed to be Anti-de Sitter (AdS) or AdS-Schwarzschild [3-5]: 

^ 2 - -{* + £ - ^ + (* + £ - ^)~V + r^ )3)) (i.i) 

where I := \J — 6/A5 is the AdS curvature radius, Mo is the black hole mass, and dtt 2 K 3 j is the maximally symmetric 

(constant curvature) 3-space with K = —1, or +1. The brane trajectory in the bulk, (i, r) = (t(r),r(r)), is 
determined by the junction condition [6]. As usual, we impose the reflection symmetry with respect to the brane. 
Then, we obtain the effective Friedmann equation on the brane as [4, 5] 

where a and p are the brane tension and energy density of the matter on the brane, respectively, and f = dr/dr 
with r being the proper time on the brane. The final term is proportional to the mass of the bulk black hole is often 
called the "dark radiation" since it behaves as the ordinary radiation. Geometrically, it comes from the projected 
Weyl tensor in the bulk, denoted commonly by E^ v [7]. If we apply Eq. (1.2) to the real universe, the values of er, I 
and M are constrained by observations of the cosmological parameters [8] . 
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When the bulk ceases to be pure AdS-Schwarzschild, or when there exists a dynamical degree of freedom other than 
the metric, the parameter M is no longer constant in general, but becomes dynamical. For instance, it is the case 
of the so-called bulk inflaton model [9-13], or when the brane radiates gravitons into the bulk [15]. In particular, in 
[10], the dynamics of a bulk scalar field is investigated in the context of the bulk inflaton model under the assumption 
that the backreaction of the scalar field on the geometry is small, and it is found that there exists an interesting 
integral expression for the projected Weyl tensor in terms of the energy-momentum tensor of the scalar field. This 
suggests the existence of a local conservation law in the bulk that directly relates the dark radiation on the brane to 
the dynamics in the bulk. 

In this paper, we investigate the case when there is non-trivial dynamics in the bulk, and clarify the relation between 
the bulk geometry and the dynamics of the brane. We focus on the case of isotropic and homogeneous branes, hence 
assume the existence of slicing by the maximally symmetric 3-space as in Eq. (1.1). In this case, we can derive a local 
energy conservation law in the bulk, in analogy with spherical symmetric spacetimes in 4-dimcnsions [16]. Then, this 
conservation law can be used to relate the brane dynamics with the geometrical properties of the bulk, especially with 
the projected Weyl tensor in the bulk. 

The paper is organized as follows. In Sec. II, we derive the local energy conservation law in the bulk and discuss the 
general property of the bulk geometry and cosmology on the brane. We show that there exists a unique generalization 
of the dark radiation that is directly related to the local mass in the bulk. We also find that there exists another 
conserved current associated with the Weyl tensor, as a non-linear version of what was found in [10]. In a vacuum 
(Ricci-flat) spacetime, the local charge for this current is found to be equivalent to the local mass. Let us call this 
the Weyl charge. The difference between the local mass and Weyl charge is given by the linear combination of 
certain components of the bulk energy-momentum tensor, and the projected Weyl tensor that appears in the effective 
Friedmann equation on the brane is indeed given by this Weyl charge. Thus we have a unique decomposition of the 
projected Weyl tensor term into the part due to the bulk mass that generalizes the dark radiation term and the part 
due to the bulk energy-momentum tensor. In Sec. Ill, as an application of the conservation law derived in Sec. II, 
we consider a simple null dust model and discuss the black hole formation in the bulk. We identify the location of 
an apparent horizon and analyze possible trajectories of the brane in the bulk. We show that the brane stays always 
outside of the apparent horizon of the black hole as long as the brane is expanding. In Sec IV, we summarize our 
work and mention future issues. 



II. LOCAL CONSERVATION LAW IN A SPACETIME WITH MAXIMALLY SYMMETRIC 3-SPACE 

In this section, we discuss the general property of a dynamical bulk spacetime with maximally symmetric 3-space, 
and consider cosmology on the brane. First, we derive a local conservation law in the bulk, as a generalization of 
the local energy conservation law in a spherically symmetric spacetime in 4-dimcnsions [16]. Namely, we show that 
a locally conserved energy flux vector exists in spite of the absence of a timelike Killing vector field. This enables us 
to define a local mass in the bulk spacetime. We also show that there exists a conserved current associated with the 
Weyl tensor. This gives rise to a locally defined Weyl charge. It is shown that the Weyl charge and the local mass 
are closely related to each other. 

Next, we introduce the brane as a boundary of the dynamical spacetime. The effective Friedmann equation is 
determined via the junction condition and it is shown that the local mass corresponds to the generalized dark radiation. 
Finally, we show that the projected Weyl tensor on the brane is uniquely related to the local mass. 

A. Local conservation law 

We assume that the bulk allows slicing by a maximally symmetric 3-space. Then, the bulk metric can written in 
the double-null form 

ds 2 = | '- dudv + r(u, v) 2 dnf K 3) , (2.1) 

where we refer to v and u as the advanced and retarded time coordinates, respectively. In Appendix A, the explicit 
components of the connection and curvature in an (n + 2)-dimensional spacetime with maximally symmetric n-space 
are listed. 

The 5-dimensional Einstein equations are given by 

Gab + A 55afc = n\T ah + S ab S(y - y ), (2.2) 
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where the indices {a, b} run from to 3, and 5, and A5 and are the 5-dimensional cosmological constant and 
gravitational constant, respectively. The brane is introduced as a singular hypersurface locates at y — yo, where y 
denotes a Gaussian normal coordinate in the direction of the extra dimension in the vicinity of the branc, and S a b 
denotes the energy-momentum tensor on the brane. The spacetime is assumed to be reflection symmetric with respect 
to the brane. 

First, we consider the Einstein equations in the bulk. They are given by 



A— [log 

- 3 - 



r 2 $ 

2f \u r \v 



log — H- + 4^- 



(*-*) 



'lij = ^iTij - r 2 7 l3 A 5 , 



(2.3) 



where 7^ is the intrinsic metric of the maximally symmetric 3-space. 

Now, we derive the local conservation law. We introduce a vector field in 5-dimcnsional spacetime as 



^ 2 



Id Id 

r „ dv r tU du 



From the form of the metric (2.1), we can readily see that £ a is conserved: 

V^?£ a ;a = (V=g?) a = 2V7((^ 3 r«),, - (r\v),u) = 0, 



(2.4) 



(2.5) 



where 7 = det jij. Note that, for an asymptotically constant curvature spacetime, the vector field £ a becomes 
asymptotically the timelike Killing vector field —(d/dt) a . 
With this vector field £ a , we define a new vector field, 



ga = ^bf b a 



where 



1 



Tab — T ab o^SSafc- 



(2.6) 



(2.7) 



Using the Einstein equations, the components of the vector field S a are given by 

4V=gS" = l\r*{K-zj\ V7, 

Z L J ti 



r z K-$) 



(2.8) 



Then, we have the local conservation law as 

5% = 0. (2.9) 
Since £ a is conserved separately, the conservation of S a implies that we have another conserved current S a defined by 



ga ._ ^o Tb a l =S a + A 5 f 



Thus we have the local conservation law for the energy-momentum tensor in the bulk. 
From Eqs. (2.8), we readily see the local mass corresponding to S a is given by [16] 



(2.10) 



M 



:= (K-$y 



(2.11) 



where the factor 3/2 in the original expression for S a is eliminated for later convenience. Alternatively, corresponding 
to S a , we have another local mass that excludes the contribution of the bulk cosmological constant, 



M :=M- -A 5 r 4 



(k - 1 



:A 5 r 4 



(2.12) 
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In what follows, we focus on the matter part M, rather than on the whole mass M. It may be noted, however, that 
this decomposition of M to the cosmological constant part and the matter part is rather arbitrary, as in the case of 
a bulk scalar field. Here we adopt this decomposition just for convenience. For example, this decomposition is more 
useful when we consider small perturbations on the static AdS-Schwarzschild bulk. 

We note that, in the case of a spherically symmetric asymptotic fiat spacetime in 4-dimensions (hence K — +1 and 
with no cosmological constant), this function M agrees with the Arnowitt-Deser-Misner (ADM) energy or the Bondi 
energy in the appropriate limits. 



B. Local mass and Weyl charge 



From the 5-dimensional Einstein equations (2.3), we can write down the local conservation equation for M in terms 
of the bulk energy-momentum tensor explicitly as 



M tV = -Kjr 3 (T u v r, u -T\r, v ) 
M,„ = |^r 3 (r„r„-TV 



or in a bit more concise form, 



dM = ^n 2 5 r 3 (T u v r, u dv + T v u r, v du - T v v dr^j . 



(2.13) 



(2.14) 



Using the above, we can immediately write down two integral expressions for M given in terms of flux crossing the 
u = constant hypersurfaces from v\ to v 2 , and flux crossing the v = constant hypersurfaces from u\ to u 2 , respectively, 
as 



of" 2 / \ 
M(v 2 ,u)-M(v 1 ,u) = -k 2 5 j dvr 3 (T u v r, u -T v v r, v ) 

2 f U2 / \ 

M(v, u 2 ) - M(v, u 1 ) = -4J dur 3 [T\r, v - T>,„J 



(2.15) 



Finally, let us consider the Weyl tensor in the bulk. In the present case of a 5-dimensional spacetime with maximally 
symmetric 3-space, there exists only one non-trivial component of the Weyl tensor, say C vu vu . The explicit expressions 
for the components of the Weyl tensor are given in Appendix A, Eqs. (A7). Using the Bianchi identities and the 
Einstein equations, we have [26] 



C. 



abed 



be: 



where 



Ja 



be 



2(n- 1) 



ln+2\T c[a . <b] t 

From this, we can show that there exists a conserved current, 

Q a = rl b n c J bca - Q a - a = 0. 
where t a and n a are a set of two hypersurface orthogonal null vectors, 



(2.16) 



(2.17) 



(2.18) 



l * = \l§(- r ,vdv) a 



n„ 



V 2 r„ \du 
i = 1/ — [ — — 



2 r, 



dv 



The non-zero components are written explicitly as 

QU = _ rJ VU v ; QV = _ rJ VU u 



(2.19) 



(2.20) 
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and we have 

(r 4 C vu vu ) ^ = r A J\ v , (2.21) 

f r 4r< vu\ _ 4 ju 

\ 1 ^ra I ' ° uu • 

\ / ,u 

These are very similar to Eqs. (2.8). It is clear that r 4 C vu vu defines a local charge associated with this conserved 
current, that is, the Weyl charge. 

Using the Einstein equations, we then find that the Weyl charge can be expressed in terms of M and the energy- 
momentum tensor as 

4 2 

r 4 C vu vu = 3M+j (QG\ - G%) = 3M + ^ r 4 (&T\ - T%) . (2.22) 

This is one of the most important results in this paper. As we shall see below, the Weyl component C vu vu is directly 
related to the projected Weyl tensor £^„, and hence this relation gives explicitly how the local mass M and the local 
value of the energy-momentum tensor affects the brane dynamics. 



C. Apparent horizons 

As in the conventional 4-dimensional gravity, the gravitational dynamics may lead to the formation of a black hole 
in the bulk. Rigorously speaking, the black hole formation can be discussed only by analyzing the global causal 
structure of a spacetime. Nevertheless, we discuss the black hole formation by studying the formation of an apparent 
horizon. 

In 4-dimensions, an apparent horizon is defined as a closed 2-sphcrc on which the expansion of an outgoing (or 
ingoing) null geodesic congruence vanishes. Here, we extend the definition to our case and define an apparent horizon 
as a 3-surface on which the expansion of a radial null geodesic congruence vanishes. Note that 'radial' here means 
simply those congruences that have only the (v, u) components, hence an apparent horizon will not be a closed surface 
if K = 0. 

The expansions of the congruence of null geodesies forming the u = constant and v = constant hypersurfaces, 
respectively, are given by [16] 

1 •„ 1 $ 1 •„ 1 $ , 

Pu = -^ a , a = — , Pv = --v a , a = -- . 2.23 

2 2r r u 2 2r r „ 

Naively, if $ = 0, one might think that both p u and p v vanish. However, from the regularity condition of the 
metric (2.1), we have 

_4^pi>0. (2.24) 

Hence, it must be that r u = or r .„ = 0, if $ = 0. If $ = r iV = 0, we have p u — and an apparent horizon for the 
outgoing null geodesies is formed, whereas if $ = r „ = 0, we have p v = and an apparent horizon for the ingoing 
null geodesies is formed. 



D. Brane cosmology 

We now consider the dynamics of a brane in a dynamical bulk with maximally symmetric 3-space [3]. The brane 
trajectory is parameterized as (v, u) — (v(t), u(t)) . Taking r to be the proper time on the brane, we have 

"'-uv = -l, (2.25) 



on the brane, where u = du/dr and so on. The unit vector tangent to the brane (i.e., the 5-velocity of the brane) is 
given by 
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and the unit normal to the brane is given by 



n 



d 



d 



dv du 



2^ \vX \v ( . 



(iidv — v duj . 



The components of the induced metric on the brane are calculated as 

dx a dx h 



<2V 



dy^ dy 



;9ab ■ 



(2.27) 



(2.28) 



where p, v run from to 3 and y 11 are the intrinsic coordinates on the brane with y° = t and y' 1 — x l (i — 1,2,3). 
Then the induced metric on the brane is given by 



-dr 2 + r(r) 2 dnf K ^, 



(2.29) 



The trajectory of the brane is determined by the junction condition under the Z 2 symmetry with respect to the 
brane. The extrinsic curvature on the brane is determined as 



where is assumed to take the form 



St = <&&g.(-p,p,p,p) - aSt ■ 



(2.30) 



(2.31) 



with a and p being the tension and energy density of the matter on the brane, respectively, as introduced previously, 
and p being the isotropic pressure of the matter on the brane. Substituting the induced metric (2.29) in Eq. (2.30), 
we obtain 



r. u u = -- 



r, v v = - 



f (p + a)+J 



2 



(2.32) 
(2.33) 



where H = r/r. Multiplying the above two equations and using the normalization condition (2.25), we then obtain 
the effective Fricdmann equation on the brane: 



M 



(2.34) 



We see that M is a natural generalization of the dark radiation in the AdS-Schwarzschild case to a dynamical bulk. 
For a dynamical bulk, M varies in time. The evolution of M is determined by Eq. (2.14), and on the brane it gives 



M = M^v + M^ u u 

2 A[ T vv (l4(p + a)-Hy-T uu (l K l(p + a)+H)u 2 



2 4 

3 5 



- -4r 4 HT\ . 



This result is consistent with [12, 15]. From the Codacci equation on the brane [7], 

D V K\ - D^K\ = K 2 5 T ab n»q\ , 



(2.35) 



(2.36) 



where is the covariant derivative with respect to q^ u and K^ v is the extrinsic curvature of the brane, we obtain 
the equation for the energy transfer of the matter on the brane to the bulk, 



p + 3H(p + p) = 2 (-T vv v 2 + T uu ii 2 ^ . 



(2.37) 



Equations (2.34), (2.35) and (2.37) determine the cosmological evolution on the brane, once the bulk geometry is 
solved. These equations will be applied to a null dust model in the next section. The case of the Einstein-scalar 
theory in the bulk is briefly discussed in Appendix B. 

Now we relate the above result with the geometrical approach developed in [7] , in particular with the E^ v term on 
the brane. The projected Weyl tensor 



(2.38) 
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has only one non-trivial component as 

E TT = C abcd n a n c v b v d = 4 C uvuv u 2 v 2 = -C vu vu . (2.39) 

Using Eq. (2.22), this can be uniquely decomposed into the part proportional to M and the part due to the projection 
of the bulk energy-momentum tensor on the brane. We find 



E T 



-^ + i(tf ( -«r.)_-^ + |(rt-«r.). (2.40) 



If we eliminate the M/r term from Eq. (2.34) by using this equation, we recover the effective Friedmann equation 
on the brane in the geometrical approach [7], 



rr comes from the projection 

case by 



where Trr comes from the projection of the bulk energy-momentum tensor on the brane and is given in the present 



T#> = - T\ . (2.42) 

Finally, from the brane point of view, it may be worthwhile to give the expressions for the effective total energy 
density and pressure on the brane. They are given by 

p(tot) _ ^(branc) _|_ ^(bulk) ^ p (tot) _ ^(brano) _|_ ^(bulk) ^ (2.43) 

where 

^^brane) = 3 [^2 („ + „)] ' , ^(brane) = l^p + ,) ( p _ ff + 2 p) , 

^(buik) = 3M ^ 4p(hulk) = M + 1 2 j_u fVu + 2f , (2 M) 

V V rj \ 1) 11 ' 



where M is given by Eq. (2.11) and f a b is defined by Eq. (2.7), and both contain the contribution from the bulk 
cosmological constant. It may be noted that, unlike the effective energy density, the effective pressure contains a part 
coming from the bulk that cannot be described by the local mass alone. The contracted Bianchi identity implies the 
conservation law for the total effective energy-momentum on the brane: 

^(bulk) + 3H ^(bulk) + p (bulk)^j = _^(bra„c) _ 3 R ^(brane) + p (bra„c)^ _ ^ ^ 

This is mathematically equivalent to Eq. (2.35). However, these two equations have different interpretations. From 
the bulk point of view, Eq. (2.35) is more relevant, which describes the energy exchange between the brane and the 
bulk, whereas a natural interpretation of Eq. (2.45) is that it describes the energy exchange between two different 
matters on the brane; the intrinsic matter on the brane and the bulk matter induced on the brane. The important 
point is, as mentioned above, that the pressure of the bulk matter has contributions not only from the local mass 
but also from a projection of the bulk energy-momentum tensor, which makes the equation of state different from 
p(buik) _ p(t>uik)^ - i e ^ j-^at of a simple dark radiation. 



III. APPLICATION TO NULL DUST MODEL 



In this section, by using the local mass derived in the previous section, we discuss the bulk geometry and brane 
cosmology in the context of a null dust model. Especially, we pay attention to the gravitational collapse due to the 
emission of energy from the brane. Namely, we consider an ingoing null dust fluid emitted from the brane [15, 17, 18]. 



A. Set-up 



The energy- momentum tensor of a null dust fluid takes the form [24] , 

Tab = Hitah + H2n a n b , 



(3.1) 
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where £ a and n a are the ingoing and outgoing null vectors, respectively, introduced in Eqs. (2.19). If we require that 
the energy-momentum conservation law is satisfied for the ingoing and outgoing null dust independently, we have 

where /(u) and g(it) are arbitrary functions of w and it, respectively, and have the dimension (G5 x mass) . We assume 
the positive energy density, i.e., f(v) > and g(u) > 0. Thus, the non-trivial components of the energy-momentum 
tensor are 

T - T - ^ 

1,17 — ^,3 1 m — 3 • v^ - ^/ 

To satisfy the local conservation law in an infinitesimal interval (u, u + du) and (v,v + dv), we find that the intensity 
functions f(v) and g(u) have to satisfy the relation, 

m& =5H(f) • (3-4) 

In general, if both f(v) and g(u) are non-zero, it seems almost impossible to find an analytic solution that satisfies 
Eq. (3.4). Hence we choose to set either f(v) = or g(u) = 0. In the following discussion, we focus on the case that 
g(u) = 0, that is, the ingoing null dust. 



B. Bulk geometry of null dust collapse 

For g(u) = 0, Eqs. (2.14) give 

M <v = \ K l — f(v), M, u = 0. (3.5) 

The second equation implies M = M{v). Substituting Eq. (3.3) into the Einstein equations (2.3), we find 

— = e F ^"\ (3.6) 
r, v 

where the function F(v) describes the freedom in the rescaling off the null coordinate v. This equation is consistent 
with Eq. (3.4). Thus, we obtain the solution as 

* = r,„e F <"> =K+^-^; M(v) = l«g jf" dv e F ^ f(v) + M , (3.7) 

where we have assumed that f(v) = for v < v n , that is, v is the epoch at which the ingoing flux is turned on. For 
dcfinitcncss, we assume that the bulk is pure AdS at v < v and set M = in what follows. 
Transforming the double-null coordinates (v,u) to the half-null coordinates (v,r) as 

r tU du = dr — r tV dv, (3.8) 

the solution is expressed as 

ds 2 = -^{r,v)e- 2F ^dv 2 + Ae- F ^dv dr + r 2 dfl 2 {K3) , (3.9) 

where $ is given by the first of Eqs. (3.7). This is an ingoing Vaidya solution with a negative cosmological constant [15, 
17]. For an arbitrary intensity function f(v), this is an exact solution for the bulk geometry. Note that if we re-scale 
v as dv — > dv — e~ F dv 1 f(v) scales as f(v) — > f(v) = e~ 2F f(v), which manifestly shows the invariance of the solution 
under this rescaling. 

An apparent horizon for the outgoing radial null congruence is located on the 3-space satisfying 

$ = r,„ = 0, while r. u = finite. (3.10) 
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This gives 




(3.11) 



The direction of the trajectory of the apparent horizon is given by 

dv 2(r 4 + M£ 2 ) 6(r 4 + MP) ' [ ' ' 

Thus, for f{y) > 0, dr/dv is positive, which implies that the trajectory of the apparent horizon is spacelike. 

For the case of K = +1 or K = 0, the apparent horizon originates from r = 0, while it originates from r — £ for 
K = — 1. A schematic view of the null dust collapse is shown in Fig. 1. We assume that the the brane emits the 
ingoing flux during a finite interval (bounded by the dashed lines in the figures) and no naked singularity is formed. 
For all the cases, the causal structures after the onset of emission are very similar. The spacelike singularity is formed 
at r = 0, but it is hidden inside the apparent horizon. 



(1) 


r > 0, 


$ > 


(2) 


r > 0, 


$ < 


(3) 


r < 0, 


$ > 


(4) 


r < 0, 


$ < 



C. Brane trajectory in the bulk 

In the null dust model, using Eq. (2.25), the proper time on the brane is related to the advanced time in the bulk 
as [18] 

4± = e , w i±^H. (313) 

To determine the appropriate sign in the above, we require that the brane trajectory is timelike, hence v > 0, and 
examine the signs of v± for all possible cases: 

-> v+ > 0, V- < 0. 
-> v+ < 0, v- < 0. 
-> v+ > 0, w_ < 0. 
-> u + > 0, w_ > 0. 

From these, we can conclude the following. For an expanding brane, r > 0, the brane exists always outside the 
horizon, <I> > 0, and v is given by v + . On the other hand, a contracting brane, r < 0, can exist either outside or inside 
of the horizon. Thus, if the brane is expanding initially, the trajectory is given by v — u+, and it stays outside the 
horizon until it starts to recollapse, if ever. If the brane universe starts to recollapse, which is possible only in the 
case K = +1, by continuity, the trajectory is still given by v — i> + , and the brane universe is eventually swallowed 
into the black hole. 

From the above result, we find 

r — \/ r 2 + $ 

r, u ii = r-r, v v= < 0. (3.14) 

Using Eq. (2.32), this gives an upper bound of the Hubble parameter on the brane as 

H< I«§( p + £r ). (3.15) 

Let us now turn to the effective Friedmann equation on the brane. For simplicity, we tune the brane tension to the 
Randall-Sundrum value, n\a — 6/£. The effective Friedmann equation on the brane is 



■> K 1 . 1 . , Mil . 
where M(r) = M(v(t)) for notational simplicity. From Eq. (2.37), the energy equation on the brane is given by 



p + 3 r -( P + P ) = -2^, (3.17) 
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r = 




= oo 



K = -\ 

FIG. 1: Causal structure of a spacetime with ingoing null dust for the cases of K = +1, and —1. In each figure, The 
(almost vertical) wavy curve represents the brane trajectory and the dotted line is the locus of the apparent horizon. The thick 
horizontal line at r — represents the spacelike curvature singularity formed there. The ingoing flux is assumed to be emitted 
during a finite interval bounded by the dashed lines. 



where /(r) = f(v(r)). From Eq. (2.35), the time derivative of M is given by 



M = -rut 



l K 2 5 (p + a)-H 



f(r)v 2 . 



(3.18) 



Thus, from Eq. (3.15), M continues to increase on the brane. 

The advanced time in the bulk is related to the proper time on the brane by v + in Eq. (3.13). Specifically, using 
the equality, 



r 2 M 



36 



(3.19) 
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on the brane, we have 

* =_ 2tU iP + <7) - H ) ■ (3 - 20) 

Note that the product fv 2 is invariant under the rescaling of v. Once /(r) is given, we can solve the system of 
equations (3.16), (3.17) and (3.18) self-consistently for a given initial condition, and determine the bulk geometry and 
the brane dynamics at the same time [15]. A quantitative analysis of the brane cosmology is left for future work. 



D. Formation of a naked singularity 

In the previous subsections, we assumed that there is no naked singularity in the bulk. However, it has been shown 
that a naked singularity can be formed in the null dust collapse [19-21]. For instance, a naked singularity exists in 
a Vaidya spacetime when the flux of radiation rises from zero sufficiently slowly. We expect the same is true in the 
present case. 

Without loss of generality, we set e F ^ — 2. We consider the following situation. For v < 0, the bulk geometry is 
purely AdS. The radiative emission from the brane begins at v = 0. We choose the intensity function as 

o \ 

f(v) = -v, (3.21) 

K 5 

where A is a positive constant. This corresponds to the self-similar Vaidya spacetime if the cosmological constant 
were absent [19]. The brane ceases to emit radiation at v = vq and the bulk becomes a static AdS-Schwarzschild for 
v > v - Thus the local mass is given by 

(v < 0) 

M(v)= I l W (0<^<«o) (3.22) 
\ Xv 2 (v < v). 

The singularity is formed at (r, v) — (0,0), and it is naked if there exists a future-directed radial null geodesic 
emanating from it. The null geodesies then form a Cauchy horizon. The trajectory of a radial null geodesic is 
determined by the equation, 



dv 2 V i 2 r 2 (v) 

Let us analyze the above equation in the vicinity of v = 0. A future-directed radial null geodesic exists if x := lim dr/dv 

v— »0 

is positive. Using L'Hopital's theorem, we obtain 

x = lim = lim — = - [K . 3.24) 

v^o v v^a dv 2 \ 3x 2 J 

It is clear that the above equation has no solution when K = or K = — 1. Hence no naked singularity is formed for 
K = or K = — 1. Therefore, we consider the case K = 1. We introduce a function, 

Q(x) =3x 3 - ^x 2 + \. (3.25) 

Then, the condition for the naked singularity formation is that Q(x) = has a solution for a positive x. The function 
Q(x) has a minimal point at x = 1/3. Therefore, the singularity is naked if 

Q(l/3) = -^ + A<0, (3.26) 



that is, 



0<A<1. (3.27) 
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Thus, the bulk has a naked singularity for small values of A, i.e., for the flux of radiation which rises slowly enough. 

Our next interest is whether the naked singularity is local or global. If it is globally naked, it may be visible on the 
brane. To examine this, we integrate Eq. (3.23). In the vicinity of v = 0, we find 



rnuii(w) =x v + + ■" ) (3.2S) 
where x is the largest positive root of Q{x) = 0; 

x = i ^1+ (l - 36 A + i 6 y/2X(l - 18 A)) 1/3 + (l - 36 A - % 6 ^/2A(1 - 18 A)) 1/3 ) (3.29) 

and 



2(5x - 1) 



(3.30) 



From the form of Q(x), we readily see that Xo monotonically decreases from 1/2 to 1/3 as A increases from to 1/18, 
and hence b is positive definite. We compare this trajectory with the trajectory of the apparent horizon. It is given 
by Eq. (3.11) with K = +1. In the vicinity of v = 0, it gives 

^M = \l^v(l-^ + ..) . (3.31) 



Since xq > ^/2A/3 for all the values of A in the range < A < 1/18, and dr^p/dv is a decreasing function of v while 
dr nun jdv is an increasing function of v, it follows that the null geodesic lies in the exterior of the apparent horizon 
and the difference in the radius at the same v increases as v increases, at least when v is small. This suggests that 
the singularity is globally naked. 

In Fig. 2, we plot the loci of the null geodesic and the apparent horizon. The result is clear. The null geodesic 
always stays outside of the apparent horizon, thus outside of the final event horizon at v = vq. Mathematically, this 
is due to the cosmological constant term in Eq. (3.23), which strongly drives the null geodesic trajectory to larger 
values of r. Thus, we conclude that the naked singularity is global and visible on the brane. The causal structure in 
this case is illustrated in Fig. 3. Investigations on the effect of the visible singularity on the brane are necessary, but 
they are left for future work. 

Finally, let us mention the strength of the naked singularity as we approach it along a radial null geodesic. Let w 
be an affine parameter of the geodesic, w = be the singularity, and the tangent vector be denoted by k a = dx a /dw. 
We examine R a bk a k b and C vu vu . From Eq. (3.3) and the Einstein equations, we have 

= y = ^ < * y - "^ . (3.32, 

r 6 \dw J r 3 \dw ) w^o x (l — xq) z 

Also, from Eq. (2.22), we have 

C vu vu = —r = —2 ► — v~ 2 cx w i-o . 3.33 

r 4 r 4 u.^o x 4 

Thus the Ricci tensor and the Weyl tensor diverge as w~ 2 and w~ T= ^o , respectively, which is a sign of a strong 
curvature singularity. 



IV. CONCLUSION 



In this paper, in the context of the RS2 type braneworld, we discussed the dynamics of the bulk and the effective 
cosmology on the brane in terms of the local conservation law that exists in the bulk spacetime with maximally 
symmetric 3-space. 

First, we formulated the local conservation law in the dynamical bulk. We found that the bulk geometry is completely 
described by the local mass M and it is directly related to the generalized dark radiation term in the effective Friedmann 
equation. We also found that there exists a conserved current associated with the Weyl tensor and the projected Weyl 
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V 



FIG. 2: The loci of the null geodesic (the solid curve) and the apparent horizon (the dotted curve) on the (v, r)-plane, scaled 
in units of the AdS radius I, in the critical case A = 1/18. Their behaviors are qualitatively the same for all the other values 
of A in the range < A < 1/18. 



tensor that appears in the geometrical approach is just the local charge for this current, and it can be expressed in 
terms of M and a certain linear combination of the components of the bulk energy-momentum tensor. 

Next, as an application of our formalism, we adopted a simple null dust model, in which the energy emitted by the 
brane is approximated by an ingoing null dust fluid, and investigated the general properties of the bulk geometry and 
the brane trajectory in the bulk. Usually, the ingoing null dust forms a black hole in the bulk. However, in the case 
oi K = +1, a naked singularity can be formed in the bulk when the flux rises from zero slower than a critical rate. 
We show that the naked singularity is global and thus it can be visible to an observer on the brane. Studies on the 
implications of a visible naked singularity on the brane is left for future work. 

Also, we found that the brane can never enter the black hole horizon as long as it is expanding. In addition, we 
found an upper bound on the Hubble expansion rate, given by the energy density of the matter on the brane, for 
arbitrary but non-negative energy flux emitted by the brane. We also presented a set of equations that completely 
determine the brane dynamics as well as the bulk geometry. 

Finally, let us briefly comment on some future issues. In this paper, we only discussed the case of null dust. However, 
this is too simplified to be realistic. As a realistic situation, it will be interesting to consider a bulk scalar field such 
as a dilaton or a moduli field. In this case, it will be necessary to solve the bulk and brane dynamics numerically in 
general. Another interesting issue will be the evaporation of a bulk black hole by the Hawking radiation and its effect 
on the brane dynamics. We plan to come back to these issues in future publications. 
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APPENDIX A: GEOMETRICAL QUANTITIES AND LOCAL CONSERVATION LAWS IN 

(n + 2)-DIMENSIONS 

In this Appendix, we give useful formulas in an (n + 2)-dimensional spacctime with constant curvature n-space, 
and generalize the expression for the local mass and Weyl charge. 
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r = 




FIG. 3: Causal structure of a spacetime with ingoing null dust when a naked singularity is formed. The wavy and almost 
vertical curve represents the brane trajectory and the dotted line is the locus of the apparent horizon. A naked singularity is 
formed at r = along the v = null line. A radial, future directed null geodesic originating from the naked singularity (the 
right-pointed thick line) stays outside of the apparent horizon and reaches the brane. 



We consider the metric in the double-null form, 

At f 

ds 2 = -^^dudv + r{u,v) 2 dn 2 (Kn) , (Al) 

where K = +1,0, or —1, corresponding to the sphere, flat space and hyperboloid, respectively. We denote the metric 
tensor of the constant curvature space as jij . The explicit expressions for the geometrical quantities in this spacetime 
are as follows. 

• Christoffel symbol 



uu 1 



ij 



log 

r $ 
27 



r $ 



-Hi 



TV. 



2r, v 



Ho 



jk jk 



(A2) 



• Riemann tensor 



7?" 
K i 



' *' ivj — 



log 



1 / $ 

-r 



2 

(log 



r <j> 

r_ u J ,u 2r u 
r $ 



(log 









)J 



7?' 











)j 




s 5 





Hi > 



15 



T>1 _ 
IX V jy 



1 ' vju ' * ujv °j 







$ 


),J 


R l jkl 


-( 



(A3) 



Ricci tensor 



Run = (lOg ) , Ryy = n^(\Og ) , 



Rii 



/, 1* .ill* ,v \ T .uv 

^* + 2( n -l)(*-*) 



rr 



1*3 



(A4) 



• scalar curvature 



• Einstein tensor 



R = — 



r-. u r.. 



( l0g -T" ) 



- 2n ; — 



+ 



n(n — 1) 



(if - $) 



(A5) 



ut ,=n(n-l) — — ^1-— j +n— , 
J r 2 $ 



G 



( l0g -T" J 



+ 2( n -l)^ 
r 



(n-2)(n - 1) 



>7ij 



(A6) 



• Weyl tensor 



n — 1 
n + 1 



log 



ijkl 



-A 



likljl ~7i 



(A7) 



From these formulas, we can show the existence of a conserved current in the same way as given in the text. 
Namely, with the timelike vector field £ a defined by Eq. (2.4), the currents S a = t; h f b a and S a = £, h T b a are separately 
conserved, and the corresponding local masses are given, respectively, by 



M 



= r n-1 



and 



(A8) 



M = M- 



-A n+2 r 



71-1 



(n - l)(n - 2) 

The w and u derivatives of M are given by the energy-momentum tensor as 

2r™ 



M v = K 2 n+2 —(T u v r,u-TV v r, v ), 
M iU = K n+2 —^-\T u r iV — T v r tU J . 



(A9) 



(AlO) 



16 



Let us now turn to the conserved current associated with the Weyl tensor. We start from the equation that results 
from the Bianchi identities [26], 



C abc / = Ja b c (All) 



where 



J abc 



2(n-l) 2 



«n+2(r c[ a;fc] + T^T 9c[ b T; a] ) • (A12) 



n n+z \ cla;Di ' {n+iy 
From this equation, we can show the existence of a locally conserved current Q a given by 

Q a = r£ b n c J bca , Q% = 0, (A13) 
where l a and n a are the null vectors defined in Eqs. (2.19). The non-zero components are explicitly written as 

We then find the following relations, 

(r n+1 C vu vu ) ^ = r n+1 J\ v , (A15) 

\ 1 *~^vu j ' " uu • 

These relations are generalization of Eqs. (2.21), and imply that the Weyl component r n+1 C vu vu is the local charge 
associated with this conserved current. 

Using the explicit form of C vu vu in Eqs. (A7) and the Einstein equations, we can relate the Weyl charge to the 
local mass. We find 

r n+i r vu _ n{n-\)M n-l +1 / \ 
T °™ " 2 n(n + l) r V G 1 2nG V 

Finally, we note that this equation implies that the linear combination of the energy-momentum tensor, 

r n+1 (T\ - 2nT"„) , (A17) 

plays the role of a local charge as well. Therefore, the behavior of this quantity is constrained non-locally by the 
integral of the flux given by the corresponding linear combination of the currents S a and Q a . Although we do not 
explore it here, this fact may be useful in an analysis of the behavior of the bulk matter. 

APPENDIX B: EINSTEIN-SCALAR THEORY IN THE BULK 

In this Appendix, we apply the local conservation law to the 5-dimcnsional Einstcin-scalar theory. We assume that 
there is no matter on the brane, but we take account of a coupling of the bulk scalar field to the brane tension. In 
this case, the energy exchange between the brane and the bulk, hence the time evolution of M, occurs through the 
coupling. 

We first consider a general bulk scalar field. Then, as a special case, we analyze the local mass on the brane for the 
exact dilatonic solution discussed by Koyama and Takahashi [13]. Finally, we clarify the relation between the local 
mass and the term which is identified as the dark radiation term in the effective 4-dimensional approach in which 
the contribution of the scalar field energy-momentum to the brane is required to take the standard 4-dimensional 
form [11]. 
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1. Set-up 



We consider a theory described by the action, 

1 „ 1 



J d 5 Xy/=l 



2 ^R-^d a 4>d a <i>-V{4>) 



WO- 



(Bl) 



For the bulk with the metric given by Eq. (2.1), the energy-momentum tensor in the bulk is given by 



T — — - 



2 ^jW^*,! 



V(<f>), 



f>,n4>,v + v(4>) 



7ij 



(B2) 



On the brane, the first derivatives of the scalar field tangent and normal to the brane are expressed, respectively, 



as 



<p := <t>, a n = -<t>,vV + 4>,uU, 

<fi := 4>. a v a = 4>. v v + <j> !U u. (B3) 
The Codacci equation (2.36) gives, via the coupling to the brane tension, the boundary condition at the brane, 

,, 1 d 



— 2^' 

In the present case, the effective Friedmann equation induced on the brane, Eq. (2.34), becomes 



*+4 



3M 



The time evolution of the local mass M on the brane is given by 



M = -l-K 2 5 r 4 H 
3 5 



1 4 4 i d 2 

36 5 Y d<f> 



(B4) 



(B5) 



(B6) 



From the brane point of view, as given by Eqs. (2.43) in the text, the effective energy density and pressure are 
composed of the brane tension and the bulk matter induced on the brane as 



(B7) 



where 



2 (T) 1 4 2 2 (T) 1 4 2 



«!p (b) 



™L K %(B)_M 1 



(B8) 



where Eq. (B4) is used. From the Bianchi identity on the brane, the conservation law for the total effective energy- 
momentum on the brane is obtained as 



}(T) 



(B9) 



The above relation is mathematically equivalent to Eq. (B6). 

As discussed after Eq. (2.45) in the text, Eq. (B9) gives the point of view from the brane, and it is naturally 
interpreted as the equation describing the energy exchange between the brane tension and the bulk matter induced 
on the brane. On the other hand, the time variation of the local mass along the brane, Eq. (B6), gives the point of 
view from the bulk, and it contains not only the energy transfer from the brane tension to the bulk (the last term) 
but also the energy flow of the bulk scalar field at the location of the brane, which is non-vanishing in general even if 
the scalar field has no coupling to the brane tension. 
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2. Dilatonic exact solution 



In the case K = 0, and for special forms of V(<p) and (j{4>), an exact cosmological solution is known, as a realization 
of the bulk inflaton model [13]. The forms of the potential and brane tension are 

K lV(4>) = ( ^ + s) A*e- 2 ^», (BIO) 



8 

kI<t(4>) = V2X e-^ bK ^, (Bll) 
where S, b and Ao are constant and are all assumed to be non-negative, and 

A = 46 2 - | . (B12) 

If S = 0, there exists a static, Minkowski brane solution [14]. In order to avoid the presence of a naked singularity, the 
dilatonic coupling b 2 is assumed to be smaller than 1/6 [13]. This implies that A is negative and it is in the range, 

2 < (-A) < | . (B13) 

The exact solution takes the form 

ds 2 = e 2W ^ (-dr 2 + e 2a (%dxW + e 2V2 ~ bK ^) dz 2^ ^ 

<j> = 0(t)+E(s), (B14) 

with the brane located at z = z and it is assumed that S(zo) = without loss of generality. The scale factor of the 
brane and the scalar field on the brane are given by 

i 

r(r) - e a(r) = (# t) ^ , e ^»*M = Hqt, (B15) 

where 



ffo = (A + -) A ^/— = 46 2 A ^/— . (B16) 

As seen from the first of Eqs. (B15), the power-law inflation is realized on the brane for b 2 < 1/6. 

Let us consider the time evolution of the energy content in this model. From the brane point of view, the time 
derivative of the brane tension p( T ) is always negative; 

* m = J^<°- < B17 » 

Thus, from Eq. (B9), for an observer on the brane, there is one-way energy transfer from the brane tension to the 
bulk matter induced on the brane. From the bulk point of view, however, the situation is slightly more complicated. 
The time derivative of the local mass (or the generalized dark radiation) on the brane, Eq. (B6), is evaluated as 



r 4 186 6 (5t 3 V3 



+ 5). (B18) 



The sign of M is determined by the sign of A/ 8 + 6. Note that the sign of A/ 8 + S determines the sign of the bulk 
potential as well, as seen from Eq. (BIO). If A/8 + 5 > 0, i.e., 5 > (-A)/8 = (b 2 /2) - 1/3, we have M > 0. Since 
M is the total bulk mass integrated up to the location of the brane, the increase in M implies an energy flow from 
the brane to the bulk. Therefore, in this case, the energy in the brane tension is transferred to the bulk scalar field 
and it flows out into the bulk. In contrast, if 5 < (— A)/8, we have M < 0. In this case, although there is still energy 
transfer from the brane tension to the bulk scalar field, the bulk energy flows onto the brane. In other words, there is 
a localization process of the bulk energy onto the brane that overwhelms the energy released from the brane tension. 
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3. Local mass and the effective 4-dimensional description 

In the bulk inflaton model with a quadratic potential [9-12], it has been shown that the bulk scalar field projected 
on the brane behaves exactly like a 4-dimensional field in the low energy limit, H 2 £ 2 -C 1, where H is the Hubble 
parameter of the brane, and the leading order correction gives the gradual energy loss from the scalar field to the 
bulk, giving rise to the dark radiation term [9, 11]. Here, we discuss the relation between the dark radiation term 
appearing in this effective 4-dimensional description and the generalized dark radiation term given by the local mass 
in the bulk. 

From the geometrical description [7], the induced Einstein equation on the brane is written as 

(A) G, V = ~4a% v + «§T$ - E^, (B19) 

where 

T$ = \ [fabq^ql - (f ab n a n b - \f ab g ah ) q„ v ] , (B20) 

is the projected tensor of the bulk energy-momentum onto the brane which includes the contribution of the cosmo- 
logical constant; see Eq. (2.7). For a homogeneous and isotropic brane, the non- vanishing components are 

rw - -f\ + 1t\, 

f^\ = \f\--f v u ~-f u v +f\. (B21) 
6 v u 

Let us decompose E^ u as 

E^ = E$+E$, (B22) 

where Ejx) is to be expressed in terms of the bulk energy-momentum tensor in such a way that the effective 4- 
dimensional description is recovered, and E$ is the part that should be identified as the dark radiation term in this 
effective 4-dimensional approach. To be in accordance with [11], we choose the components of E$ as 

-4 6) = -\4(^T\ + |T\) + ^4(T\ - 3T\) = -£<»\, (B23) 
and identify the remaining part with the dark radiation term, X, 

-E { t f =X = -E^ i i . (B24) 
In the effective 4-dimcnsional description, the Einstein equation on the brane takes the form, 

MGp V = KlTW-E$, (B25) 

where T^f ^ is the effective energy-momentum tensor on the brane, 

«2r w = - V + e$> , (B26) 

and k\ is the 4-dimensional gravitational constant that should be appropriately defined to agree with the conventional 
4-dimensional Einstein gravity in the low energy limit. In the present case of homogeneous and isotropic cosmology, 
the only non-trivial components are the effective energy density and pressure, which are given explicitly by 

nlp^ = -«Jf <«*>% = \n\f\ + \k\T\ - \4(p v u + |f\), 

4p (cS) = Ulf^\ = ~4v 2 + Uf\ + - U(-f\ + ?f\). (B27) 

6 12 4 12 8 \u u / 

The effective Fricdmann equation on the brane is written as 

K 



H +^2 



= Y 2 <° 2 \4f\ + - A 4f\ - '-4 (Hf \ + lf\) + X. (B28) 
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Applying this to a bulk scalar field with the action (Bl), we find p( eS ^ and p( eff ) are given by those of a 4-dimensional 
scalar field ip with the potential, 



2 KV> 12 yv> 16 



(B29) 



where ip — (f>. From the contracted Bianchi identity we obtain 



i(eff) 



+ 3H 



>ff) +p (eff) 



1 1 



(B30) 



Unfortunately as we can see from Eqs. (B27), there is no simplification in the energy equation in terms of the 
5-dimensional energy-momentum tensor. 

From the effective 4-dimensional point of view, what happens is the conversion of the scalar field energy on the brane 
to the dark radiation via the coupling to the brane tension. From the bulk point of view, a natural interpretation is 
to regard the local mass M on the brane as the generalized dark radiation. These two different identifications of the 
dark radiation term on the brane coincide only when the bulk is in vacuum and M is constant. Comparison of the 
above decomposition of E TT with Eq. (2.40), we find the difference between the dark radiation in the 4-dimensional 
description and the generalized dark radiation in terms of the local mass M as 



3M 



(B31) 
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